home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / IMTQLV.z / IMTQLV
Text File  |  1996-03-14  |  4KB  |  67 lines
  1.  
  2.  
  3.  
  4. ____IIIIMMMMTTTTQQQQLLLLVVVV((((3333FFFF))))                                                        ____IIIIMMMMTTTTQQQQLLLLVVVV((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      IMTQLV, SIMTQLV  -  EISPACK routine.  This subroutine is a variant of
  10.      IMTQL1.
  11.  
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  14.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee  iiiimmmmttttqqqqllllvvvv((((nnnn,,,, dddd,,,, eeee,,,, eeee2222,,,, wwww,,,, iiiinnnndddd,,,, iiiieeeerrrrrrrr,,,, rrrrvvvv1111))))
  15.           iiiinnnntttteeeeggggeeeerrrr          nnnn,,,, iiiieeeerrrrrrrr,,,, iiiinnnndddd((((nnnn))))
  16.           ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn dddd((((nnnn)))),,,, eeee((((nnnn)))),,,, eeee2222((((nnnn)))),,,, wwww((((nnnn)))),,,, rrrrvvvv1111((((nnnn))))
  17.  
  18.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ssssiiiimmmmttttqqqqllllvvvv((((nnnn,,,, dddd,,,, eeee,,,, eeee2222,,,, wwww,,,, iiiinnnndddd,,,, iiiieeeerrrrrrrr,,,, rrrrvvvv1111))))
  19.           iiiinnnntttteeeeggggeeeerrrr          nnnn,,,, iiiieeeerrrrrrrr,,,, iiiinnnndddd((((nnnn))))
  20.           rrrreeeeaaaallll             dddd((((nnnn)))),,,, eeee((((nnnn)))),,,, eeee2222((((nnnn)))),,,, wwww((((nnnn)))),,,, rrrrvvvv1111((((nnnn))))
  21.  
  22.  
  23.  
  24. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  25.      This subroutine finds the eigenvalues of a SYMMETRIC TRIDIAGONAL matrix
  26.      by the implicit QL method and associates with them their corresponding
  27.      submatrix indices.
  28.  
  29.      On INPUT
  30.  
  31.      NNNN is the order of the matrix.
  32.  
  33.      DDDD contains the diagonal elements of the input matrix.
  34.  
  35.      EEEE contains the subdiagonal elements of the input matrix in its last N-1
  36.      positions.  E(1) is arbitrary.
  37.  
  38.      EEEE2222 contains the squares of the corresponding elements of E. E2(1) is
  39.      arbitrary.  On OUTPUT
  40.  
  41.      DDDD and E are unaltered. Elements of E2, corresponding to elements of E
  42.      regarded as negligible, have been replaced by zero causing the matrix to
  43.      split into a direct sum of submatrices.  E2(1) is also set to zero.
  44.  
  45.      WWWW contains the eigenvalues in ascending order.  If an error exit is made,
  46.      the eigenvalues are correct and ordered for indices 1,2,...IERR-1, but
  47.      may not be the smallest eigenvalues.
  48.  
  49.      IIIINNNNDDDD contains the submatrix indices associated with the corresponding
  50.      eigenvalues in W -- 1 for eigenvalues belonging to the first submatrix
  51.      from the top, 2 for those belonging to the second submatrix, etc.
  52.  
  53.      IIIIEEEERRRRRRRR is set to ZERO       for normal return, J          if the J-th
  54.      eigenvalue has not been
  55.         determined after 30 iterations.
  56.  
  57.      RRRRVVVV1111 is a temporary storage array. Calls PYTHAG(A,B) for sqrt(A**2 +
  58.      B**2).  Questions and comments should be directed to B. S. Garbow,
  59.      APPLIED MATHEMATICS DIVISION, ARGONNE NATIONAL LABORATORY
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.